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In a wide variety of natural and laboratory magnetized plasmas, filaments appear as a result of 
interchange instability. These convective structures substantially enhance transport in the direction 
perpendicular to the magnetic field. According to filament models, their propagation may follow 
different regimes depending on the parallel closure of charge conservation. This is of paramount 
importance in magnetic fusion plasmas, as high collisionality in the scrape-off layer may trigger 
a regime transition leading to strongly enhanced perpendicular particle fluxes. This work reports 
for the first time on an experimental verification of this process, linking enhanced transport with 
a regime transition as predicted by models. Based on these results, a novel scaling for global 
perpendicular particle transport in reactor relevant tokamaks such as ASDEX-Upgrade and JET is 
found, leading to important implications for next generation fusion devices. 

PACS numbers: 52.35.Py, 52.35.Ra, 52.55.Fa, 94.30.cq 


Interchange instabilities are an ubiquitous phe¬ 
nomenon in magnetized plasmas where a density 
gradient coexists with a parallel force density of the 
same sign [lj. This instability leads to the formation 
of elongated structures known as filaments, flux tubes 
or prominences featuring a dipole electric polarization 
in the perpendicular plane. The resulting E x B drift 
propels the filament, greatly enhancing convective 
perpendicular transport. Examples can be found in 
disparate contexts, including astrophysical plasmas such 
as accretion discs [2] and planetary magnetospheres 
[3], where the destabilizing force is typically of the 
centrifugal type, and laboratory plasmas [4, 5], where 
the force density typically comes from magnetic pres¬ 
sure gradients. This issue is particularly relevant for 
magnetically confined fusion plasmas, as it determines 
the propagation of filamentary structures, which have 
become recognized as the dominant radial transport 
mechanism in the region between the closed magnetic 
field lines and the wall, known as Scrape-off Layer 
(SOL) [6-9]. Therefore, filamentary transport strongly 
influences the parallel/perpendicular ratio of the particle 
and heat fluxes onto plasma facing components, thus 
determining the durability of plasma facing components 
and the sputtering of impurities from the main wall. 

Basic models for filaments in fusion literature (see, 
e.g. Krasheninnikov [5]) describe how elongated struc¬ 


tures propagate as the result of equilibration of plasma 
polarization caused by an effective gravity force: taking 
a reduced MHD approach on a magnetized stationary 
plasma under some arbitrary force density F, the charge 
conservation equation can be expressed as 

V 'dt(^ V±<; 4 j = B b ' VXr + V!|J|1 ’ (1) 

where n, <f>, B, J and rrii stand for density, potential, 
magnetic field, current density and ion mass, d/dt = 
d/dt+(hVcj))/B-X and b = B/B. Polarization (lhs of Eq. 
1) is thus the result of an equilibrium between the drive 
(first term on rhs) and the parallel closure term (second 
term on rhs) accounting for the current parallel to B. In 
the SOL context, F = 2 nrriiC 2 s IR e r , stands for the effect 
of curvature and X/B (c s is the sound speed and e r indi¬ 
cates the radial direction). This can be easily generalized 
to many other forcing mechanisms such as the centrifugal 
force in Keplerian systems, gravity, etc. Equation 1 yields 
different solutions depending on the parallel closure used 
for Jy: the primary hypothesis in the SOL is that the 
filament extends along the field line to the solid wall, 
and J|| is thus limited by the sheath [7]. In this regime, 
known as “Sheath Limited” (SL), the polarization term 
is neglected and V||J|| = j^nc s e(l - exp [-e(0 - (j>f)/T e ]) 
is canceled by the drive term [7]. Here, (f>f , e, T e and L\\ 
are the floating potential, electron charge, electron tem¬ 
perature and parallel connection length. By solving Eq. 
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1, the perpendicular velocity of the filament, Vb, is then 
found to scale as 


Vb ~ 2Cs " 6*' (2) 

where p s and R are the ion gyroradius and tokamak ma¬ 
jor radius, and 5b is the filament perpendicular size. How¬ 
ever, several phenomena, such as high SOL collisional- 
ity, large X-point magnetic shear, electromagnetic effects, 
etc. may invalidate the wall connection hypothesis, effec¬ 
tively disconnecting the filament from the wall [7]. In this 
case, V|| J\\ -*• 0 and polarization cancels the drive term, 
leading to a different regime, known as “Inertial Regime” 
(IN), in which [ 1] 


Vb = c s yjp e ^ oc 5l 12 , (3) 

where p e is the pressure in the filament normalized to 
the background value. These basic models have been 
refined by the inclusion of realistic SOL elements such 
as hot ions [11] and large fluctuation amplitudes (“full-f” 
approximation) [12], leading to improved versions of Eqs. 
2 and 3: 


n 5 b 
n + n R 

(4) 

where = Ti/T e , and n is the background density. Note 
that the isothermal limit has been assumed to approxi¬ 
mate p ^ Tn. These scaling laws should be considered 
as upper boundaries, as any deviation from the pure in¬ 
terchange cross-phase between electric field and pressure 
fluctuations reduces the velocity of the blob with respect 
to predictions [13, 14], Although these models have been 
successfully compared with experiment in basic plasmas 
[15] and extensive characterization work has been made 
in tokamaks [L6-18], a direct measurement of the transi¬ 
tion between the two regimes remains to be achieved in a 
fusion relevant plasma. Nevertheless, such transition has 
been invoked to explain the formation of the SOL density 
shoulder, i.e. a substantial increase of the far SOL den¬ 
sity radial e-folding length, A„ ^ (Vr-n/n ) -1 , observed in 
many tokamaks when a certain density is exceeded dur¬ 
ing L-mode operation [19-21]. Recent experiments in 
ASDEX Upgrade (AUG) [22] linked the increase in \ n 
to an increase of the filament size and associated perpen¬ 
dicular particle transport in the outer midplane T ± fn, 
which can be up to 40% of the total transport after the 
transition according to comparison to EMC3-EIRENE 
simulations [23]. This transition takes place as the col- 
lisionality increases in the SOL at the onset of divertor 
detachment, suggesting that the aforementioned discon¬ 
nection might play a key role. Myra et al. [24] predicted 
this process in a two region model of the SOL using the 
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as the control parameter, where v e i is the electron-ion 
collision rate and fstands for the gyrofrequency of 
ions/electrons. The disconnection takes place for A > 1, 
when the characteristic parallel transport time is longer 
than the inverse of ion-electron collision frequency. 
Subsequent work indicates that this results in enhanced 
perpendicular transport as the consequence of increased 
filament velocities and creation rates [25]. Results 
from Ref. [22] showed that the transition in AUG 
coincides with A ^ 1. The relation between filament 
regime transition and the shoulder formation has also 
been shown in JET [26], where independent measure¬ 
ments of collisionality in the divertor and midplane 
regions suggested that only the collisionality in the di¬ 
vertor region could account for the disconnection process. 


The question of whether the collisionality in the di¬ 
vertor or in the midplane, Aor A mi d, are determining 
the shoulder formation is of great practical importance, 
as next generation tokamaks are foreseen to operate 
with partially detached, locally collisional divertors, 
while remaining hot and collisionless at the midplane. If 
shoulder formation is dominated by a filament transition 
induced by divertor collisionality, this phenomenon will 
not be reduced by a larger machine size, as it would ac¬ 
cording to other proposed mechanisms such as main wall 
recycling [27]. In such case, a shoulder will probably be 
formed [26], which would impact substantially main wall 
particle fluxes and erosion. To solve this conundrum, 
an experiment was designed in AUG to separate the 
contributions of A m id and Ain L-mode discharges, 
T e at the outer midplane far SOL does not depend 
strongly on the injected power Pin, and remains nearly 
constant at T emi d ^ 20 eV for a wide range of densities. 
In contrast, divertor detachment depends strongly on 
Pin , occurring at lower midplane densities for lower 
heating power. Here, A m id and A dm are defined as in 
Eq. 5, using respectively n and T e Langmuir probe 
measurements from the midplane and target plates. In 
the case of A m ,d, L\\ is the connection length from the 
measurement point to the wall, L\\ ^ irRqgs. In the 
case of Adiv, it has been estimated as the connection 
length between the target and the X-point height at 
the flux surface of the measurements L\\ ^ ^TtRqgb. 
This assumption is based on reciprocating probe and 
spectroscopic measurements of the cold region extension 
in front of the divertor at the time of the transition 
[ 12, 28]. Therefore, by realizing a series of density 
ramps with different heating powers (0, 300 and 600 kW 
of ECH power, plus 500 kW of Ohmic heating in all 
cases), A m id is kept constant while A di V covers a range 
of more than two orders of magnitude, crossing the 
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Figure 1: Collisionality parameter A at midplane and 
divertor. Colors indicate different heating powers. The 
onset of divertor detachment is indicated by solid 
vertical lines. 


Figure 2: Size of filaments vs. collisionality. Colors 
indicate Pin as in Fig. 1. Thin solid lines provided as 
guide to the eye to mark the two regimes. 


A div - 1 critical point at different values of edge density, 
n e . This is shown in Fig 1, where the collisionalities 
at the divertor and midplane for three Pin cases are 
displayed. While A m id - 0.5 for most of the n e range, 
an exponential growth of Adi V can be observed at the 
onset of the divertor detachment (highlighted by colored 
dashed lines). The onset of detachment is marked by the 
DoD > 1 [29]. Details on the diagnostic setup and data 
analysis are analogous to those described in Ref. [22]. 

Two different regimes appear approximately at each 
side of Adiv = 1, as can be seen in Fig. 2. For low 
collisionality, the size of filaments remains around one 
cm, and increases logaritmically with A div However, for 
dominant collisionality (Adiv > 1), the size of filaments 
increases by up to an order of magnitude. The data set 
shows no dependence on Pin, hence the transition is 
not determined by a n e threshold, as can be deduced 
from the different Adi V (PiN,n e ) paths seen in Fig. 1. 
Furthermore, the velocity of small filaments follows the 
SL scaling 1/6% and larger ones the IN scaling \/6/,. 
This is shown in Fig. 3, where the measured average 
filament velocity v ± and size 6b are normalized with c s 
and p s respectively as in Ref. [I ]. Also, the upper 
bounds in the two regimes are represented by plotting 
the respective expressions of Eq. 4. In order to do this, 
the local magnetic field (Bsol = 1-4 T), curvature radius 
(Rsol - 2.15 m), n e and T e values at the probe position 
have been considered. To estimate h, the isothermal 
hypothesis has been extended to probe measurements, 
and the ion saturation fluctuation level ji jSa t has been 
used as a proxy. Comparison of Fig. 3 with Fig. 2 
reveals how both groups of data points correspond to 
those featuring Adiv < 1 and A div > !'■ The horizontal 


dotted line crossing the elbow in Fig. 2, is represented 
vertically in Fig. 3, where it separates the two scalings. 
Again, the same behaviour is seen for different Pin 
values, indicating the scaling is independent of this 
parameter (and thus also of n e ). For both regimes, cold 
and hot ions have been considered (t^ = 3, as proposed 
in Ref. [11])- Filaments in the IN regime are better 
represented by the cold model. This is consistent with 
both theoretical analysis [18] and far SOL measurements 
carried out with field analyzers [30], which suggest that 
filaments cool down after the transition. Data points 
in the lower left corner of Fig. 3, corresponding to 
the smallest values of 6b, follow the IN scaling rather 
than the SL. This effect has been predicted [31] for 
SL filaments with 6b « <5* - 15 p s , which would revert 
to an IN-like scaling. Finally, the electromagnetic 
scaling proposed in Ref. [11], in which electromagnetic 
effects dominate collisionality in the V|| J|| term, repre¬ 
sents a worse fitting to the data than the SL and IN ones. 

Finally, global perpendicular particle transport is 
found to be linked to divertor collisionality. In Fig. 4, the 
density profiles remain almost constant until A div - 1 is 
surpassed. For A div > 1, the SOL width rapidly increases 
by a factor of 3, indicating a substantial increase in F ± . 
As before, all parameter scans in AUG match onto one 
curve when A di V is used as ordering parameter. Here, A„ 
is derived from Li-beam spectroscopy in the first 25 mm 
outside the separatrix, as explained in Ref. [22]. These 
results are independent on how detachment is achieved: 
If nitrogen is puffed in the divertor, local radiation cools 
down the region inducing detachment at a lower range 
of densities for a given Pin■ In Fig. 1, seeding would 
cause a shift of the dashed line to the left compared to 
the unseeded discharge at the same Pin- However, the 
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Figure 3: Filament scaling. Color code as in Fig. 1. 
Black dotted line indicates the regime separation, as 
observed in Fig. 2. Red/Black lines indicate IN/SL 
regimes, as expressed in Eq. 4. Solid/dashed indicates 
Ti - 0/ri = 3. Blue solid line indicates the 
electromagnetic regime described in Ref. [11]. 



Figure 4: Density e-folding length scaling with divertor 
collisionality. Black, blue and red circles represent Pin 
as in Fig. 1. Purple circles correspond to seeded 
discharges (empty/solid mean before/after N 2 puffing). 
Light blue stars are data points from JET. Straight 
lines mark the two regimes as in Fig. 2. 


A„ values from these discharges (empty/solid purple dots 
indicating before/after the puffing) scale exactly as the 
unseeded ones, despite the different h e and Pin range 
(see Fig. 4). Measurements in JET yield very similar 
results: X n and Aare calculated for the discharges 
presented in [26] , using equivalent lithium beam and tar¬ 
get probes data. The resulting data, displayed in Fig. 
4 as light blue stars, display a remarkable similarity to 
the ones measured at AUG. As can be deduced from Fig. 
1 this scaling is only obtained with Adi V , and not with 
Amid■ This is consistent with previous studies showing 
that local collisionality at the midplane does not strongly 
influence midplane transport [27]. Beyond magnetic fu¬ 
sion, these results can be regarded as the validation of 
a generic model for the propagation of structures result¬ 
ing from interchange instabilities in magnetized plasmas. 
In this sense, the behavior of SOL filaments is strikingly 
similar to that of filamentary structures observed in the 
Jovian and Saturnian magnetospheres by the Galileo and 
Cassini spacecraft [32, 33]: The radial density profile of 
the equatorial iogenic plasma torus in the Jovian mag¬ 
netosphere is dominated by the convection of “isolated 
interchanging flux tubes” driven by an isomorphical cen¬ 
trifugal force-driven interchange instability [34]. In this 
case, the parallel closure condition is determined by the 
conductivity at the end of the field lines, represented by 
the ionospheric conductivity at high planetary latitude 
instead of the plasma-wall interaction in the divertor re¬ 
gion. As in the SOL case, when the coupling to the 
ionosphere is lost due to low conductivity, disconnected 
filamentary structures increase their radial speed, and, 
in qualitative agreement with Galileo measurements, the 


radial gradient is substantially decreased with respect to 
regions in which the equatorial plane and the poles are 
connected [35]. 

Summarizing, the role of A di V as a control parameter 
for the transition of both filament and global perpen¬ 
dicular transport regimes has been demonstrated, after 
separating its effect from other plasma parameters such 
as n e , A m id or Pin- The velocity scaling of filaments 
changes from 1/6% to \/6f } when the critical value 
A di V = 1 is surpassed. These results represent strong 
evidence of filament propagation being governed by the 
mechanism presented in [24], in which the interruption 
of the parallel circuit due to collisionality switches 
filaments from the SL to the IN regime. This change 
in SOL turbulence is accompanied by a substantial 
increase of the perpendicular particle flux T ± in the 
far SOL, measured in two different tokamaks. The 
remarkable similarity of AUG and JET results supports 
the prediction of a shoulder formation in ITER. This has 
important implications for main wall particle fluxes and 
erosion in future fusion devices. Finally, the similarity 
between filaments driven by isomorphic instabilities in 
such disparate contexts as the SOL of fusion devices and 
planetary magnetospheres points towards a universal 
feature of transport in magnetized plasmas. 
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